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Abstract 

We study the existence and stability of static kink-hke configurations of a 5D scalar field, with 
Dirichlet boundary conditions, along the extra dimension of a warped braneworld. In the presence 
of gravity such configurations fail to stabilize the size of the extra dimension, leading us to consider 
additional scalar fields with the role of stabilization. We numerically identify multiple nontrivial 
solutions for a given 5D action, made possible by the nonlinear nature of the background equations, 
which we find is enhanced in the presence of gravity. Finally, we take a first step towards addressing 
the question of the stability of such configurations by deriving the full perturbative equations for 
the gravitationally coupled multi-field system. 
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I. INTRODUCTION 



The possibility of extra spatial dimensions [l, 2|, hidden from our current experiments 
and observations through compactification or warping, has opened up a wealth of options for 



particle physics model bui 
cosmological problems 117 



dmg 



16| and has allowed entirely new approaches for addressing 



32| . In many implementations, Standard Model (SM) fields can 



be confined to a submanifold, or brane, while in others they populate the entire extra- 
dimensional space. Common to both approaches, however, is the inclusion of bulk fields 
beyond pure gravity, either because they are demanded by a more complete theory, such 
as string theory, or because they are necessary to stabilize the extra-dimensional manifold. 
Thus, a complete understanding of the predictions and allowed phenomenology of extra 
dimension models necessarily includes a comprehensive consideration of the configurations 
of these bulk fields, the simplest of which are real scalars. Indeed, 4D Poincare invariance 
allows for these new bulk fields to acquire nontrivial static configurations along the extra 
dimensions. 

Static one-dimensional scalar configurations with a node (where the field vanishes) are 
known to localize wave functions of other fields near that node. In the context of extra 
dimensions, these kink-like scalar backgrounds can be used for example to localize bulk 



fermions near either boundary 



33|-|36| , allowing for interesting constructions of flavor models. 



They can also affect the localization of other scalar or vector fields leading to a field theoretic 



description of fat branes (see for example the constructions in 37H39[|). Kink-like scalar 
configurations are a particularly interesting case to consider because the boundary conditions 
make it possible to obtain non-trivial general results regarding both the existence and the 
stability of such configurations, at least in the case of one fiat extra dimension without 
gravity {40! . 41]. In this paper we build on these previous results and extend them as far 
as possible to the case with a gravitating (warped) extra dimension. In the presence of 
gravity the kink-like configuration cannot fix and stabilize the interbrane distance {4^. It 
is therefore necessary to assume the existence of at least one additional stabilizing field, 
coupling either directly or gravitationally to the kink field. We will opt for the latter and 
introduce additional non-interacting scalar fields. At least one of these additional fields 
must be given a monotonic profile in order to stabilize the size of the extra dimension (i.e., 
interbrane modulus) 43 1. 
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The plan of the paper is as follows. In section II we review the results for kink-like 
backgrounds in a flat extra dimension with no gravity. We then generalize these to include 
a warped gravitational background, but with no gravitational backreaction from the kink- 
scalar itself. In section IV we consider the coupled system of multiple scalar fields in the 
presence of gravity, and as a special case consider a kink field with Dirichlet boundary 
conditions and another scalar whose purpose is to stabilize the whole configuration. We 
write the background equations for such a system and show graphically how non-linearities 
allow a given action to have multiple static solutions. 

In the final section we take the first steps towards studying the question of the stability 
of these static configurations by deriving the complete set of equations for scalar and grav- 
itational perturbations around a given static background. The general procedure is quite 
complex and involves extended theorems of oscillation theory appropriate to the type of 
eigenvalue problem we are lead to, namely a matrix Sturm-Liouville problem. We therefore 
reserve a complete study of the stability of the general system for future work. 



II. KINKED SCALARS IN FLAT EXTRA DIMENSIONS 



In 



40 



4l| a 5D flat scenario including one real scalar field with an arbitrary scalar 
potential was studied and the general conditions for the existence and perturbative stability 
of static, nontrivial, background scalar field configurations were presented. In this section 
we briefly review the main results and slightly extend the discussion of the energy densities 
of different kink configurations. 

Consider a real scalar field in 5 dimensions (labeled by indices M, N, . . . = 0,1,2, 3, 5) 
with a flat background metric, and defined by the action 

'1 



S= (Fx 



(1) 



The extra dimension is compactified on an S1/Z2 orbifold with the scalar field (l){x,y) being 
odd under Z2 reflections along the extra coordinate x^^y (i.e. (f){x,y) = —4>{x, —y))- Here 
the orbifold interval is defined as [0, ttR], with its size ttR assumed to be fixed. The potential 
V{(f)) must then be invariant under the discrete symmetry (p — )■ —0, and is chosen to have 
at least two degenerate global minima at = ±f, with v^O. To simplify notation, we will 
also choose the potential to vanish at = 0. 
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FIG. 1: Profiles in tlie extra dimension interval [0,7ri?] of different static configurations of the 
Dirichlet scalar field (p, defined by the scalar potential V{(j)) = + |A|(/>^ (/i^ = 2M^, 

A = lil/^^, 7ri?=8.6375M~^). The solutions with nodes in the interval (dashed curves) are unstable, 
while the stability of the nodeless and trivial solutions (solid curves) depend on the parameters of 
the model. 



Under these conditions, it was shown in 



40| that there will always be static solutions. 



nontrivial along the extra coordinate satisfying the (static) field equation 

, 



dV 

90 



(2) 



where a prime denotes a derivative with respect to y. The profiles of these solutions, satisfy- 
ing Dirichlet boundary conditions, resemble that of a kink solution patched to an anti-kink 
in the middle of the interval. The possible solutions were classified in two groups, namely 
those with nodes in the interval (multiple kink-antikink solutions patched together) and 
those with no nodes, vanishing only at the end-points of the orbifold (see Fig. [T]). It was 
shown that all static kink solutions with nodes are perturbatively unstable, whereas the sta- 
bility of nodeless solutions depends on the parameters of the model in a particularly simple 
way. 

The Dirichlet solutions of Eq. (j2]) with no nodes in the interval form a continuous one- 
parameter family of functions. A simple choice for the parameter is the amplitude A of the 
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FIG. 2: Nodeless static configurations of tlie kink scalar field (/>, defined by the scalar potential 
V{(j)) = —\\^'^\(jP' + |A|(/>'* (/i^ = 2M^, A = \M^^). Configurations with different amplitudes are 
solutions to different physical problems, corresponding to different stabilization radii of the extra 
dimension. The vertical dashed line indicates the minimal radius Rc^ below which nodeless solutions 
do not exist with this potential. 

solution, i.e., the maximum value of the nontrivial solution (pAiv)- Solutions with different 
amplitudes A generally vanish at different points along the extra dimension, which corre- 
spond to different possible orbifold radii R (see Fig. [2]). However, in order to obtain the 
stability condition for these solutions it is extremely useful to consider the full family of 
solutions. 

The value of the 4D effective energy density of a given static solution 0a (y) is 

E{A) = (^^<p'/ + F(0a)) dy , (3) 

where T{A) is the length of the solution in the extra dimension. This can be conveniently 
rewritten as an integral over using properties of Eq. ([2]) and its solutions 0a (l/) 

h ^V(4,) - V(A) ^ 



We are now equipped to state the general results of 40|, |4l| in a slightly modified, although 
more revealing, version: 



Proposition 1 A static solution to equation with S > nodes inside the orbifold in- 
terval is always unstable. 



Proposition 2 A static, nodeless solution (pA^in) to equation with amplitude A^,, and 
associated energy density E{A^:) is stable if 



dE 
dA 



< . (5) 

A=A, 



This is a powerful result since it means that given any scalar potential V{(j)) we immediately 
know which of the nontrivial nodeless solutions (pA will be stable or unstable, without the 
need to actually know explicitly their analytic form. 

With this result it is possible to understand the vacuum structure of any single scalar field 
theory with Dirichlet boundary conditions in 5D when the metric along the extra dimension 
is fiat. Possible static solutions consist of the trivial solution (0) = (which may or may not 
be stable), kink- like solutions with nodes in the interval (which are always unstable), and 
kink-like solutions without nodes in the interval (some stable and some unstable, depending 



on condition (15])). As remarked in 



40 



4l| . the trivial solution may be the true vacuum 
solution even in the case of a negative mass term — in the 5D potential, as long as 

the inequality |yU^| < |1/-R^| is preserved. Therefore, for a given orbifold radius R, many 
different perturbatively stable vacuum solutions are possible, and it is necessary to identify 
which one is the true vacuum of the theory. 

The true vacuum of the theory will depend on the size of the radius R. This can be 
seen as follows: Without loss of generality, one may define the energy density of the trivial 
solution to be zero by choosing the 5D potential V{(j)) to vanish at = 0. It was shown 



m 



40 



41| that there is a critical radius Rc below which nontrivial nodeless solutions do 



not exist (see Fig. [2]). The energy density associated with the critical nontrivial nodeless 
solution will be either positive or exactly zero, so that the transition from one vacuum to 
another can be either second order or first order, as one varies the radius R. 

In Fig. |3] we show an example of a simple setup defined by the scalar field potential 
V{(f)) = _ i|A|04 + 11^106^ ^j^i^ ^2 _ 4^2^ ^ _ ^^-i ^ ^ O.QM-\ In 

the right panel, the energy density of two static solutions is plotted as a function of R, 
showing clearly that below a critical radius -Ri only the trivial solution is possible and above 
a critical radius R2 only the kink solution is possible. For Ri < R < R2, both solutions 
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FIG. 3: Profiles (left panel) in the extra dimension interval [0, vri?] of the two possible stable static 
configurations of the Dirichlet scalar field (p, defined by the scalar potential V{(j)) = — ^|/i^|(/>^ — 
\\X\(t)^ + (with /i2 = 4M2, a = 4M-1, ^ = 0.6M-^ and vri? = 1.368M-1). In the right 

panel, we show the energy of the two stable solutions as a function of vri?, and it is seen how the 
absolute stability of coexisting static configurations is determined by the size of the radius R. The 
dots show critical points where the scalar perturbations contain a massless mode. 

are perturbatively stable. At the radius i?* the two solutions are degenerate, marking the 
transition from one true vacuum to another {(ptm for R < R* and 0kink for R > R*)- From 
this we see that the inverse length scale 1/R plays the role of an order parameter of a phase 
transition, much like temperature T in finite temperature field theory. For a very small 
radius R (analogous to high T) the system is stable only around its trivial solution, with all 
symmetries restored. As the radius increases (analogous to T decreasing) the system can 
undergo a phase transition, which could be of either first or second order. The analogy with 
temperature, however, is not meant to be taken literally. For whereas the temperature in 
any 4D effective cosmology must be monotonically decreasing for most of its history, the 
orbifold radius R could in principle increase, decrease or oscillate on very long time scales, 
depending on the dynamics of the stabilization mechanism (which we have so far ignored). 

III. KINKS ON A WARPED BACKGROUND 



We now extend previous investigations to the case of a scalar field in a warped extra 
dimension, while neglecting any backreaction on the warping from the scalar field itself. 
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In this case one includes the effects of the curved metric along the extra dimension on 
the scalar field solutions while still ignoring the dynamics of the gravitational sector. We 
therefore consider the action 

"1 



S 



(TXy 



-9 



2' 



{dM4>) dN<P - V{(t>) 



where the form of the metric is now taken to be 



(6) 



(7) 



and where a{y) is the warp-factor and 7^1, the 4D metric on slices of constant y. The purpose 
of considering scalar field configurations on a fixed background is to explore whether our 
previous results continue to hold in the presence of a warped background in a regime where 
we still have semi-analytical control over the solutions. We postpone a discussion of the 
full dynamical problem, including the backreaction on the metric due to the presence of the 
scalar field, until the next section. 



A. Kink Scalar in an AdS^ Background 

In the original Randall-Sundrum (RS) model {q], the metric takes the form ([7]) with 
a{y) = k\y\ and 7^1, = r]^^, where k has dimensions of mass and is related to the 5D 
cosmological constant of AdS^. In this background any static nontrivial field configurations 
are solutions of 

dV 



0" - 4k(j)' 



d(j) 



. 



Scalar perturbations around this kink background, ip{x,y) = 
posed as 

^{x,y) = j2^^:'\M^\y) 



(8) 

(x, y) — (l){y), can be decom- 

(9) 



such that the normal modes (p^\x) and ^p^\y) are solutions of 

(^^□y^^ + m^^, = 



where 



2 ^ dfV 



and *^^-*n = rj^'^d^dy. Taking the derivative of the kink equation 



(10) 

(11) 

gives 
(12) 



where we have defined ^pu = 4>'- Thus ipM is a massless solution (m^ = 0) of the perturbation 
equation f|TT]) . ahhough it satisfies mixed boundary conditions rather than the Dirichlet 
boundary conditions imposed on (py. 

At this point we are aheady able to state a new result of this work, which is an extension 
of the previous result related to the impossibility of having stable kink solutions with nodes 
inside the interval. Suppose that happens to have 6 nodes inside the interval. We 
have just shown that 0' = i^m will solve the equation for a massless excitation, but with 
mixed boundary conditions. Since has S nodes, 0' = ipM must have 6 + 1 nodes inside the 
interval. The following inequalities relating the eigenvalues Aj? for the Dirichlet case and the 



eigenvalues for a general mixed boundary condition case 4^ hold from Sturm-Liouville 
theory 

Xn < < ■ (13) 

Since we have A^^^ = (i.e. the eigenvalue of the solution with 6 + 1 nodes), we can 
immediately deduce that the mass-squared of the lowest excitation of the Dirichlet problem 
must be negative since A^^^ < A^^^ = with 6 > 1. 

Proposition 3 In a warped background on a slice of AdS^, any static solution to equa- 
tion with 6 > nodes inside the interval is always unstable. 

However, for nodeless static solutions (when 6 = 0) the results for the fiat case obtained 



m 



3, 



4l| ] cannot be extended here. Lacking a general stability condition, we will instead 
propose a weaker sufficient stability condition for these and other more generic solutions in 
the next subsection. 



B. Kink Scalar on a General Background 



In a general warped background with metric ansatz ([7]) the equation for a static scalar 
background configuration is 

- . (14) 



In this situation, although we have been unable to extend the stability theorems found 
earlier, we are still able to find a general sufficient condition for perturbative stability of the 
background configurations. 
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Small perturbations around the background (f){y) may be defined as in ([9]). The spectrum 
of these perturbations consists of solutions to the eigenvalue problem 



'^yiVi) = Vy{y2) = (15) 

v>y - 4aV; - [f^\y) - ^y^^y^] ^y = o. (16) 

A useful form of this equation is obtained by performing a change of variables e^^^^dy = dz 
and defining a{z) = -|ln(J(2;)) and W{z) = ix^{z)e~'^"^y''^^'^ to yield 

^-y-(W^(-)-O=0. (17) 



45| . For any function f{z) 



To proceed, we make use of the following integral inequality 
such that /(a) = f{h) = 0, and with n nodes within the interval [a,b], there exists p E TZ 
such that 

/ e-P/"W//W dz>{n + l)ev/p? . (18) 

J a 



Applied to f[T7j) . this implies 



f ^-P[^+w{z)\ > + l)e^ , (19) 

J a 



the logarithm of which yields 

ml>- ln[(n + l)e^^ --\n( e^P^^^"^^'^^ d^ (20) 

which is a lower bound for the eigenvalues in terms of the background quantities cr{z) and 
fi{z) (which are contained in J and W). In the case of the lowest eigenvalue we have 

ml>- He^) - - In r e-f^^+^^'^^ dz) (21) 

P P \J a J 

and so a sufficient condition for perturbative stability (m^ > 0) is 

^-p[Jl+w{z)] < _ (22) 



We may formulate this explicitly in terms of the warp factor a{z) so that finally, a static 
solution (p^z) of (JH]), obeying Dirichlet boundary conditions, is stable if 

J a 
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where = ^X- _ and where the actual value of p is that which extremizes the right 

</.{2) 

hand side. This is a sufficient condition, but not a necessary one. In order to demonstrate 
how effective this weaker stability condition can be, we now turn to a simple example in 
which the condition can actually be evaluated. 

Consider a flat metric, where a{y) = 0, and a trivial background scalar configuration, 
i.e., (j){y) = 0, but where the 5D scalar potential is allowed to have a tachyonic mass. In this 
case equation f[T6l) becomes 

- {P' -ml)^y = 0- (24) 
If (fy has Dirichlet boundary conditions, the solutions to this problem are 



= sin(^m2 -/i2 y) (25) 

where — p^ = r^-n'^ jh'^ and L = 6 — a is the size of the extra dimension. The mass of 
the lightest mode is = fi'^ + tc'^/L'^ and so the condition for stability is mo < , which 
means that the bulk scalar mass /i^ can be negative, but not arbitrarily so: 

/i' > -vrV^' • (26) 

Therefore in this case (where a' = a" = 0), our sufficient condition (|23l) becomes 

e-P"' dz< , (27) 

J a 

which leads to 

yu'>-ln(— ^). (28) 



The value of p that extremizes this bound is p = vrL^/e, and so our weaker bound is 

> —J, ■ (29) 

This result is a factor of 27r/e weaker than the exact bound (!26|) . Nevertheless this result 
is nontrivial as it clearly demonstrates that it is possible to have negative bulk masses and 
retain a stable system.^ 



^ The stability conditions of the trivial vacuum in the presence of negative bulk mass terms in an extra- 
dimensional scalar field theory have been analyzed and generalized to general warped backgrounds in 



46| 
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IV. KINKS IN GRAVITATING WARPED EXTRA DIMENSIONS 



So far we have examined static scalar field configurations in a fixed background. We have 
found that some of the results that were shown to hold in a fiat extra-dimensional background 
continue to hold in a fixed warped background, and we have found useful generalizations 
of other results. We now want to include the dynamics of the gravitational sector and 
explore how these results can be extended when the gravitational backreaction is included. 
Therefore we now seek nontrivial static field configurations in which the warp factor has its 
own dynamics determined by the 5D Einstein equations. 

As soon as we include a dynamical gravitational sector, we are required to worry about 
stabilization of the extra dimension. In the above discussion we assumed that the extra 
dimension was stabilized and that the dynamics of the stabilization mechanism were frozen 
out. Here we want to include the backreaction of any matter fields on the 5D metric, and 
so we must include the dynamics of stabilization. A natural question to ask is whether 
the kink fields of interest could provide a stabilization mechanism. Unfortunately, in 42| it 
was shown that when one considers static solutions for both the warp factor and a single 
scalar field, the lightest scalar perturbative mode (the radion) will be tachyonic whenever 
the derivative of the scalar profile vanishes inside the interval. In other words, the system 
is unstable whenever the scalar field profile passes through an extremum in the bulk. This 
means that if we insist on obtaining a nontrivial configuration for a single scalar field with 
Dirichlet boundary conditions, we are guaranteed to obtain a tachyonic radion and the 
extra dimension will be unstable. To address this issue we will add extra scalar fields whose 
purpose will be to stabilize the radion as in [43 1. 



The resulting system becomes considerably more difficult to analyze than the case with 
only one bulk scalar field, particularly with regard to questions about stability. On the other 
hand, the case with three or more scalar fields is formally no more difficult to analyze than 
the case with only two scalar fields. Hence we will keep our treatment general to include an 
arbitrary number of scalar fields Xa (a, 6 = 1, . . . , A/"), although when we consider particular 
examples below, we will specialize to the case with only two scalar fields (a kink field and 
a non-kink field). For simplicity we will assume throughout that the scalar fields are only 
coupled gravitationally. 
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We therefore consider the 5D action for gravity and TV free scalar fields 

M3 



S 



2 



^ / d^Xy/^ [7^ - 2A] 



£ ^^'''^(^^^^'^)(^^^«) - ^(^") - E ^^(Xa)S{y - yd , (30) 

_a=l i=l,2 

where M^, = (SttG)"^/^, G is the 5D Newton's constant, TZ is the 5D Ricci scalar, and A is 
the 5D cosmological constant. The full scalar potential in the bulk is W{xa), and the brane 
potentials are Xi{Xa)- As before, we take the 5D line element of the form 



ds^ = e-^''^y^-f^^{x)dxf'dx'' - dy^ 



(31) 



where 7^1, is the induced metric on the 4D hypersurfaces of constant y, which foliate the 
extra dimension. The 5D Einstein and field equations are 



A 



2M3 

A (4)7^ 



1(^1 1 2 

* \ a i=l 



Ai(Xa)5(y - Vi 



a 



dW 

dxa 



dX 



i=l,2 



(32) 

(33) 
(34) 



where '■^^T^ is the 4D Ricci scalar associated with the induced 4D metric 7^j,, which we have 
left arbitrary. The boundary conditions for the system are determined by Israel junction 
conditions at each brane. These are obtained by integrating the equations of motion over 
an infinitesimally small interval across each brane, giving 

1 



3M3 
dX 



X^iXa)\ 



l^'a]y^ = li™, IXaiVi + e) - x'aiVi ' ^)] = JT^ 
" €^0 dXa 



(35) 
(36) 



Vi 



These yield A/" conditions on each brane, which is exactly the number of data that need to 
be specified in order for equations and (IMI) to form a well-posed problem. 

Note that the above boundary value problem consists of a system of coupled nonlinear 
differential equations. Finding solutions analytically for such a setup is highly unlikely, 
although it is still possible to proceed in the opposite direction, i.e. given a particular 
analytical solution one can obtain the setup from which it originates. To do so, one relies 
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on the powerful method of the superpotential 47H49| . which can be useful even for two or 



more scalar fields (see, for example, 50| in the context of soft-wall models). However, even 
if one solution is constructed in this way, there is no guarantee that this is the only solution 
with the same action. We will now describe how to look for all possible solutions of a given 
action using a combination of numerical and graphical techniques. 

A. Multiple Solutions 

Whenever there is more than one static solution to the above boundary value problem with 
the same action, we say that multiple solutions exist. In general, the bulk scalar fields can 
have Dirichlet boundary conditions, Neumann boundary conditions or more general mixed 
boundary conditions. Here we focus on the case where we have one kink field (obeying 
Dirichlet boundary conditions), with the remaining A/"— 1 fields Xa having Neumann or mixed 
boundary conditions. When the profiles of these extra fields are monotonic, they will tend 
to stabilize the extra dimension, whereas if their profiles have vanishing derivatives inside 

42I. Despite this subtlety, we 



the interval, they will tend to destabilize the extra dimension 
will generically refer to the non-kink fields as "stabilization" fields. 

To find solutions we proceed as follows: we specify the Lagrangian in the bulk and on one 
of the branes, and we numerically solve an initial value problem to determine the profiles of 
the fields along the extra dimension. Dirichlet boundary conditions are imposed on the kink 
field at the initial brane by demanding that it vanish there. For this to hold, we assume 
the kink field has a sufficiently heavy brane mass so that it decouples from the stabilization 
fields on the branes. As a result, the kink field disappears from the junction conditions 
f l35|) -f l36|) . which then yield only A/" conditions on the initial brane. This leaves A/" + 1 
initial conditions that need to be specified, which we take to be the boundary values for the 
derivatives 0', x'l • • • xV-i' ■ After solving the initial value problem for a given choice 

of initial conditions, we impose Dirichlet boundary conditions on at the final boundary by 
locating the second brane at a point where the profile of vanishes. In general, the profile 
will vanish at several points along the extra dimension, and one may study kinks with the 
desired number of nodes by choosing the location of the second brane accordingly. Here, 
as in the flat case, we are primarily interested in nodeless kink solutions, and we therefore 
place the second brane at the first zero of the profile function. 
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We now have a solution to a boundary value problem whose boundary conditions on the 
second brane are not yet known. We parameterize the brane potential on the second brane 
X2{Xa) in terms of P parameters aj, (for example, the brane tension E2, the brane mass term 
ml of each scalar, the quartic coupling of each scalar, etc.) 

A2(Xa) = /(S2, ml, ml, m^, ...) . (37) 

Then the junction conditions (I35l)-(l36l) at the second brane {i = 2) give M linear equations 
for the P unknowns ab- By evaluating the fields on the second brane, and using the param- 
eterization in ( 137|) . we then invert the Af junction conditions to determine the ab. If this is 
possible, then the solution to our initial value problem is also a solution to a corresponding 
boundary value problem. From this we see that we must have P > Af in order to guarantee 
that the field configuration we obtained is the solution to a corresponding boundary value 
problem. If P = Af, the ab are uniquely determined, and there is a unique Lagrangian for 
which the above field configuration is a solution. On the other hand if P > Af, some of the 
ab are arbitrary and so there is a family of solutions for these final-boundary conditions. In 
that case there is a family of Lagrangians which yield the obtained field configuration, and 
one can proceed by focusing on one member of this family. If P < Af, the linear system of 
parameters may be overdetermined, in which case the obtained field configuration is not 
a solution to any corresponding boundary value problem. 

We can find additional solutions by changing the initial-boundary conditions and repeat- 
ing the above process. Note that by freely varying the field derivatives (0', Xi - ■ ■ xV-i) 
the initial brane and determining the remaining quantities from the junction conditions, it 
is possible to leave the initial-brane potential unchanged. This is necessary in order that the 
action remains unchanged (it is not sufficient because part of the action is determined by the 
final-brane potential). A solution and the resulting final-boundary conditions (the ab) are 
then found as before. Since each set of initial shooting values yields a set of ab, each ab is a 
function of the Af initial-boundary derivatives. Each ab therefore defines an A/'-dimensional 
surface whose level-surfaces can be projected onto the (p' {yij-x'Svi) parameter space (which 
is an A^-dimensional space). In the above construction there are P such quantities, and so P 
level-surfaces intersect at every point in this parameter space, representing one solution for 
this action. The question of whether multiple solutions exist for the same action is equiva- 
lent to the question of whether the same P surfaces simultaneously intersect at more than 



15 



one point in the parameter space. 

We will now show how this works in two simple examples. In both cases, we will consider 
a kink field in addition to just one stabilization field x, with no interaction terms among 
them in the scalar potential. In both examples there will be regions of parameter space in 
which two distinct static configurations are possible for the same action. 



B. Example 1: Quartic Potential 



In both of the following examples we consider a Lagrangian for two scalar fields 

1=1,2 

+ I 9''''{dMX)dNX - U{x) - ^^(x)S{y - Vi) , (38) 



i=l,2 



where (j) is the kink field and % is the stabilization field with potentials 

U{x) = \mlx^ (39) 
HX) = Q^.^x' + S.) . (40) 

The fact that the second brane potential for x is parameterized in terms of two parameters, 
fil and S2, will allow us to find unique solutions to the boundary conditions on the second 
brane. The junction conditions ( |35|) - (!36|) become 

^'im) = (-1)^"'^ (^l^^h^y^) + s.) (41) 
x'im) = {-ly-'l^^h{y^) ■ (42) 

On the second brane {i = 2) these can be inverted to give 

S2 = -6M^a{y2) + x'{y2)x{y2) (44) 

so that once we determine the fields on the second brane, we can extract the boundary 
conditions (and therefore Lagrangian) to which those fields are a solution. 

The only things left to specify are the bulk potential for the kink field and the initial- 
boundary conditions. In this first example, we take the kink potential to be 

V(0) = -^mj02 + iA0^ (45) 
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Taking the initial brane to be located at y = 0, we find solutions to the initial-boundary 
value problem at this brane with Dirichlet boundary conditions imposed on the field (j). 
Examples of nodeless solutions to the initial value problem are shown in Fig. |H To ensure 
that (j) obeys Dirichlet boundary conditions on the second brane, we locate the second brane 
at the first point (other than y = 0) where the profile of vanishes (the vertical dashed 
lines in Fig. H]). ^ Once the position of the second brane is identified, the final-boundary 
conditions are determined from fH3|) and (jS]). By varying the initial shooting conditions, 
(f)[ = 4>'{yi) and x'l = x'{.yi)i ^^id repeating this process of finding solutions, identifying the 
location of the second brane, and determining the final-boundary conditions, we generate 
level-curves of n\ and S2. These are plotted in Fig. [5l Notice that most /i^ contours cross 
each E2 just once, signifying that there is a single solution for the corresponding action with 
the kink potential of Eq. fH5|) . However, some contours cross each other more than once 
(see, for example, the circles in Fig. [51 ) Furthermore there is only a finite region in the 
parameter space where solutions exist. If either or are increased sufficiently, the 
solution to the initial value problem blows up. In that case the boundary value problem has 
no solution, since a second boundary where = does not exist. It is therefore possible to 
scan the entire allowed space and examine whether multiple solutions with the same 

action exist. 



C. Example 2: Higher-Order Potential 

In this second example we take a slightly more complicated kink potential 

^(0) = -\mlct^' - + • (46) 

The other potentials and boundary conditions are the same as in the previous example, the 
only difference being the dynamical evolution of the system due to the new potential V{(j)). 
We choose this potential because, contrary to the potential in our first example, in the limit 
of weak gravity and fiat spacetime, it leads to multiple solutions to the same boundary value 



problem 40|, |41| due to the nonlinear nature of the equations. 



^ For certain initial conditions, the profile of 4> will blow up before it vanishes for a second time. When this 
happens the initial conditions used do not lead to a solution of our boundary value problem. 
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FIG. 4: Profiles of the scalar backgrounds (/'(y) and x(y) as well as the warp factor o"(y), showing 
the two possible solutions (panels A and B) to the same boundary value problem defined by the 
physical parameters m\ = -O.SM^, fij = -0.25M^, fi^ = -8M^, Si = -2M^, Sa = 0.52Mf, 
ml = 0.5M2, A = 2M~\ and A = 0. 

In our more general setting, including gravity and a stabilization field, we find numerically 
that there exists more than one solution for the same Lagrangian in a large portion of the 
parameters space. In Fig. E] we show two such solutions. Note that these solutions would 
be extremely difficult to discover by randomly guessing initial-boundary conditions. To be 
more methodical we follow the same procedure as before to find level-curves of the final- 
boundary conditions, shown in Fig. [71 Again, solutions to a particular action will be given 
by the intersection of the appropriate contours for the brane mass squared and brane 
tension S2. As can be seen, there are regions in which some contours intersect at more than 
one point, showing that multiple solutions for the same action are possible as expected. In 
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S2 - - (^)7^ 




FIG. 5: Level-curves of ^2 ^2 in the parameter space for example 1 with = — 0.5M^, 

= -0.25M2, El = -2M^, = 0.5M^, A = 2A4"^ and A = 0. Circled are two points in the 
0'i"Xi parameter space with the same values of /i^ and S2, corresponding to two solutions with the 
same Lagrangian (plotted in Fig. U]) . 

particular, we again circle two such points, corresponding to the solutions plotted in Figs. Ei 
As an interesting remark, note that both of these particular solutions happen to lie near the 
region of parameter space where the 4D cosmological constant vanishes. 



V. STABILITY OF SOLUTIONS 



Having shown how different nontrivial static field configurations exist in warped extra 
dimensions, the next question to ask is whether these solutions are stable. As we reviewed in 
section [Til the case of flat extra dimensions there exist 4l| techniques for determining the 
stability of such solutions. Indeed, for certain potentials, in that case perturbative stability 
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FIG. 6: Profiles of the scalar backgrounds 0(y) and xiu) well as the warp factor showing 
the two possible solutions (panels A and B) to the same boundary value problem defined by the 
physical parameters m\ = -O.SM^, = -0.25M2, = -5M^, Si = -2M^, Sa = 0.56Mf, 



0.5M^, X = 2M-\ i = 6M-4, and A = 0. 



can be determined analytically. Unfortunately, in the case of warped extra dimensions, the 
question of stabihty is complicated by the presence of multiple scalar fields and their coupled 
dynamics. Here we begin to study the perturbative stability of these kinked configurations. 
We derive the linearized equations and reformulate the problem in terms of a matrix Sturm- 
Liouville problem. However, the full analysis requires matrix Sturm-Liouville methods which 
we omit and leave for future work. 

We begin by expanding the metric to first-order. Instead of the coordinates in ([7]), in this 
section it will be more convenient to choose coordinates so that the metric takes the form 



(47) 
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FIG. 7: Level-curves of ^2 ^2 in the (j)'-x' parameter space for example 2 with = — 0.5M^, 
fij = -0.25M2, Si = -2M^^ = O.SM^, A = 2M-^, ^ = 6M~4 and A = 0. Circled are two 
points in the (p'l-x'i parameter space with the same values of ^2 ^2, corresponding to two 
solutions with the same Lagrangian. These solutions are plotted in Fig. [6l 

Working in the generalized longitudinal gauge (see appendix lA] for details), we introduce 
scalar perturbations $ and \1/ and write the perturbed metric as 

ds^ = a\y) [(1 + 2$(x, y))-f^,{x)dx^dx'' - (1 + 2^(x, y))dy^] . (48) 

Next we expand the A/" scalar fields to first-order in small perturbations ^a{x,y) 



Xa{x,y) = Xa{y) +^a{x,y) 



(49) 



and compute the linearized Einstein equations, yielding A/" + 1 dynamical equations for 
A/" + 1 scalar fields {Af fundamental scalars and one graviscalar, or radion). Since only Af 
of these equations are independent, the Einstein constraint equations are used to eliminate 
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one of the scalar fields in terms of the others (see Appendix |X] for details). The resulting M 
independent equations are 



' 2 xV ^XN 



4a 



3M3 



2 



a=l 



/ dW 










X 







(4) 



dW 



dXa 



(50) 
(51) 



where "H = ^, (^^D = '^^^d^dy, and in equation (EI]), as throughout, we have assumed 
that there are no direct couplings between the 5D scalar fields in the scalar potential (in 
Appendix El we derive the general form of these equations when couplings between the fields 
are included). These dynamical equations can be written more compactly as 



um + Vim = vl^ 



(52) 
(53) 



where ^ has suppressed discrete indices which run over the M —1 fundamental scalar fields, ^ 
is the graviscalar, the are y-dependent differential operators (i.e., linear differential oper- 
ators having y-dependent coefficients and acting only on functions of y) also with suppressed 
discrete indices, and □ = 'y^'^d^di, is the 4D wave operator. 

The boundary conditions are determined by integrating the equations of motion across 
each brane. Integrating equations fISU]) and fIFIl) . these are found to be 



1 dW 



X'm dxM 



m' 



4a^n 



1 dW 



x'm dxM 



dxl 



Aa_ 

Me 



2 



a=l 



f dW 




Xa 


dW 




\dXa 


X 


x!n 


dxM 



2 a' 



dW 



Vi 



dXa 



(54) 



(55) 



These boundary conditions can be put in the form 



i\x,yi) = Bi{yi)m{x,yi) + B2{yi)i{x,yi) 



(56) 
(57) 
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where A12 and -81,2 are functions of defined via (155|) . and we have used (156|) in (15^ to 
obtain fl57|) . 

The plan now is to perform a separation of variables in order to obtain a Sturm-Liouville 
eigenvalue problem, and then to analyye this eigenvalue problem to determine stability of 
the system. Because the 5D equations of motion of the scalar perturbations are coupled, 
the correct separation of variables ansatz is a coupled one 

^W(a;,y) = ¥^'\y) (58) 
&\x,y) = ^^''Xy) u^^\x) , (59) 

where u^"'\x) is the ra*'* 4D Kaluza-Klein physical mode and '$^\y) and C,y"\y) are the wave 
functions. Plugging this ansatz into equations ( 152|) and ( 153|1 leads to the following coupled 
equations 

^y{y)Duix)+u{x)Vf^yiy) = u{x)V%{y) (60) 
UvWx) + <x)'^lUy) = uixm^i'yiy) . (61) 
The separation of variables thus yields a 4D wave equation for u{x) 

nu{x) + mlu{x) = (62) 
and a system of two coupled differential equations 

Vf^yiy) - ml^yiy) = VlUv) (63) 

with boundary conditions for the profiles 

%{y,) = A,{y,)^y{yi) + A2{yi)Uyi) (65) 
e;(y.) = B,iy,)^y{y,) + B2{yi)Uy^) ■ (66) 

The system of equations fl63|) and flM|) constitute an eigenvalue problem. The stability of the 
static background around which we have added scalar perturbations therefore depends on 
the existence, or absence, of a negative eigenvalue associated with a solution to eqs. ( 163|1 
and (IMl). 

This situation is somewhat unusual, since generally the Kaluza-Klein eigenvalue problem 
arising from dimensional reduction consists of a single second order differential equation, 
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which can be put in standard Sturm-Liouville form. Analyzing that Sturm-Liouville eigen- 
value problem is straightforward, since in particular it is known that the eigenvalues are 
bounded from below, and that the eigenfunction corresponding to the smallest eigenvalue 
has no zeros within the interval. Therefore, the question of stability in practical terms be- 
comes the search for a solution to the Kaluza-Klein equation such that it contains no nodes. 
Its associated eigenvalue will be the lightest possible eigenvalue and, if positive, the system 
will have no classical instabilities. 

In the present case, however, the Kaluza-Klein problem is a system of coupled differential 
equations. Consequently, matrix Sturm-Liouville techniques are required. In order to ana- 
lyze stability further, one must extend the theory of oscillations and the concept of nodes 
of solutions to a higher dimensional problem. Such an analysis, although rather involved, is 
underway, and will be presented in a future work. 



VI. DISCUSSION AND OUTLOOK 



Braneworld theories generally lead to scalar degrees of freedom that propagate in the 
extra-dimensional bulk. Understanding the vacuum structure of these models in the presence 
of bulk scalar fields is therefore a prerequisite to fully appreciating their phenomenological 
possibilities. Furthermore, bulk scalars may provide a useful way to localize fermions and 
build braneworld models purely with field theory (e.g., fat branes and soft walls). 

In this work, we have studied the vacuum structure of braneworld models with one 
warped extra dimension and multiple bulk scalar fields. In particular we have focused on 
static configurations along the extra space coordinate where one of the fields-with Dirichlet 
boundary conditions-acquires a nontrivial kink-like profile. To find these solutions one needs 
to solve both the Einstein and the scalar field equations. In the limit of a fiat 5D metric 
and weak gravity such solutions are known to exist, and the problem of finding all possible 
static confi gura tions as well as determining their perturbative stability has been addressed 



and solved 
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4l| . Here we have built upon this previous work to determine how warping 



along the extra dimension effects the existence and stability of these kink-like solutions. 

When considering a fixed warped background, it was sufficient to look for nontrivial 
solutions for a single scalar field. In this case, neglecting any backreaction of the scalar field 
on the gravitational dynamics, we found that such kink-like solutions do indeed exist. As 
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in the case of a flat extra dimension, we were able to prove that any kink-hke solution with 
nodes in the bulk is unstable. Thus we have focused on nodeless kink solutions and the 
trivial solution. However, in contrast to the flat case, in the presence of warping we were 
only able to find a sufficient condition for dctcrming the stability of these solutions. We 
were therfore unable to analytically determine stability for nontrivial solutions in a warped 
background, even when that background is fixed (e.g., in the Randall- Sundrum model with 
no backreaction) . Instead we were forced to determine stability numerically. 

Including the dynamics of the gravitational sector forces the inclusion of additional scalar 
fields whose purpose is to stabilize the size of the extra dimension. In that case we were 
again able to find nontrivial kink-like configurations, except now for a coupled multiple-field 
system. We have described a general graphical technique to find all possible static config- 
urations of the background equations with one kink scalar field and an arbitrary number 
of additional "stabilization" fields. The technique amounts to generating solution surfaces 
by varying the shooting parameters needed to solve the coupled system of equations. This 
technique also allows us to look for multiple solutions with the same action. We have demon- 
strated how to implement this technique in two simple examples, where we considered one 
kink field and one stabilization field in the presence of gravity. As in the fiat case, when 
the potential for the kink field is a higher-order polynomial (leading to higher-order non- 
linearity in the field equations), we found that multiple solutions may exist for the same 
action. Interestingly, however, we also found multiple solutions for the same action when 
the kink potential was a fourth-order polynomial, which differs from the result obtained in 
a flat background. 

We have addressed the issue of stability only partially. We have derived the full 5D 
perturbative equations, including gravitational perturbations, for multiple scalar fields in 
the presence of a warped extra dimension. The system of equations constitute a matrix 
eigenvalue problem, which must be analyzed using an extension of the usual theorems coming 
from oscillation theory or Sturm-Liouville eigenvalue problems. Such techniques exist in the 
mathematical literature but due to the complexity of the task, we have left the numerical 
analysis of the general case for a later work. 
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Appendix A: Scalar Perturbations in the Generalized Longitudinal Gauge 

In this appendix we derive the linearized 5D Einstein and field equations for scalar per- 
turbations in the bulk. We linearize around a background metric of the form 

ds'^ = a^{y) {-ff,^{x)dx^dx'' - dy'^) . (Al) 

The background Einstein and field equations in these coordinates are 

-H' - ^a^ = (j2 Ix'a' + Wixa) + ^ X.{Xa)Ky - y^)\ (A2) 



6 \^2^'' 

\ a 



- ^a'+^^ = ^\y^ '^Xl - a^W{xa) 1 (A3) 



6 12 

~dXa dXa 



x: + mx'a - -a'Y. ^^^y - = , (A4) 
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where 1-i = ^, "^^^T^ is the 4D Ricci scalar with respect to the background 4D metric 7^j,. 
To first-order in scalar perturbations, the 5D metric can be written 



+ 2Bix, y\^dx^dy - (1 + 2^{x, y))dy^] 



(A5) 



where | indicates a covariant derivative with respect to the 4D slices of the bulk. Choosing 
to work in the generalized longitudinal gauge, we set S = -E = 0, and the linearized metric 
simplifies to 



ds^ = a^{y) [(1 + 2$(x, y))-i^y{x)dx^'dx'' - (1 + 2^(x, y))dy'^] . 
We also expand the scalar fields to first-order 



(A6) 



(A7) 



where the fields Xa obey the background equations of motion (lA2p - (lA4p above and C,a{x,y) 
are small perturbations. The linearized Einstein and field equations are 



2$ + ^ = 



(A8) 
(A9) 



(^)n(2<l> + ^) - 4$" + + 8^2* + 4H(^' - 3$') + -(^)7^<l' 

3 
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* a \ 



dW 
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* a \ 



dW 



dxadxb 
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dXa 

2x'> - - 4$' + 6-Hvl>) 



(AlO) 



(All) 



(A12) 



where (^^D = Y'^d^d^ is the 4D wave operator. Applying the constraint equation f lA8p to 
equations (1A9l) -( IAT2l) . and making use of the background equations (lA2l) -( IA4l) . yields 



(A13) 
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(A14) 

^al (A15) 



dXa 



^. (A16) 



We obtain a 5D wave-like equation for \E' by subtracting (]A14p from (lAlSP to give 



(4)q^ _ ^" _ 97^^' _ ( 47^' + 127^2 _ _(4)^ ] ^ 



3M3 



9ty 



(A17) 



Equations (]A17p and ( 1A16P comprise JV+I dynamical equations for the A/'+ 1 perturbation 
variables \I' and ^a- However, since these variables are connected through the constraint 
flAlSp . only A/" of them are independent. Therefore we may use flA13P to eliminate one of 
the variables in terms of the others. Choosing to eliminate the TVth scalar field, ^^v") in terms 
of ^a<j\f and equation ( lAlSI) gives 



1 f^^' 



(A18) 



(Note that we cannot eliminate \l/ in terms of the scalar fields ^a, since this requires an 
integration over unknown functions. This is understandable since doing so would amount to 
reducing the problem to one in flat spacetime, which ought to be impossible.) Substituting 
this into ( 1A17I) and (1A16I) and rearranging gives 



(4) 



1 dW 



x'm ^xm 



2 xV ^XX 



A 2 A'-l 



3M3 



E 

a=l 



/ dW 




X'a dW 




\dXa 


X 


XV ^XM 





(A19) 



(4) 



b=l 



X'b 



-3 Xa 



2xV dxMdXa 



dXbdXa 



XV dxMdXa 



dW 



dXa 



+ 



2xV ^XAfdXa 



^ . (A20) 
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We may write this more compactly as 



N-l 



(^)OT+Pfi[/-5^(P|Ua =0 

a=l 

wne. + 5^(i)IU6-raa^ =0, 



(A21) 
(A22) 



b=i 



where 

ma 



m - 2a^ — 



1 dW 



4a2 / dW 



3M3 \ dXa 



1 9„ - I 12H^ + AH' - l^^^n - Aa^n^ 
-J \ 2 xV 



1 dW 



3)ab 



-Sabid^ + Sndy) + Mab 



dxbdxa 



+ 



X'b d^W 



XV dxAfdXa 



ma = -3 xi 



2xV dXNdXa 



dW 



dXa 



+ 



2xV dx^dxa 



(A23) 



If we suppress the discrete indices in equations flA21l) and flA22l) they take an even simpler 
form 



(A24) 
(A25) 



This is the generic form that the scalar perturbation equations of motion take for A/" — 1 
coupled scalar fields and a graviscalar. 
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